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0\ ; Abstract 

I We investigate high density state of SU(2) QCD by using Lattice QCD simulation 

T-H ■ with Wilson fermions. The ratio of fermion determinants is evaluated at each step 

. of the Metropolis link update by Woodbury formula. At /? = 0.7, and k = 0.150, we 

I calculate the baryon number density, the Polyakov lines, and the energy density of 

' gluon sector with chemical potential ^=0 to 0.8 on the 4^ x 12 lattice. Behavior of 

^~~| . the meson propagators and diquark propagators with finite chemical potential are 

y—i , also investigated. 

o ■ 

: 

^ : 1 INTRODUCTION 

^ ' High density state of the strongly interacting matter is attracting much attention [|l| and 
^ ■ one of the main targets of this workshop. Numerical simulation based on the Lattice 
^ . QCD is the established method to investigate the high temperature state of the strongly 
interacting matter; However, because of the well known problem that chemical potential 
makes action complex, the progress of the lattice QCD in the finite density has been rather 
slow. Indeed, after the first dynamical quark simulation with the chemical potential was 
done for SU(2)[^, to our knowledge, only few full SU(3) QCD calculations had been tried. 
Recently, Fodor and KatzQ proposed the nobel method to draw the critical line on the 
T — /i plane and to find the point where first order phase transition turns to the crossover. 
Their method is based on the Lee- Yang Zero, therefore, it seems to work well only on 
the critical line and it seems still out of scope to investigate physics across the phase 
transition. 

Quark in the real world is SU(3) fundamental representation; SU(2)-QCD and quark 
in the adjoint representation are simple toys for the theorist. However, due to the re- 
cent progress in analytical investigations, we can hopefully obtain some information on 
real SU(3)-QCD through the investigation of the finite density region of the "QCD-like" 
theories 0. The QCD-like theories, such as SU(2)-QCD, quark model in the adjoint 
representation and QCD at finite isospin density, are expected to have less difficulties in 
numerical analyses. In these years, there are indeed high activities in Monte Carlo calcu- 
lations with dynamical quark of such kinds of models H]. Furthermore, recent analyses 
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on the color superconductivity suggest the possible realization of SU(2) part of the color 
SU(3) as a residual interaction in the color superconducting state f^. In this paper, we 
report our recent work on the SU(2)-QCD finite density states with Wilson fermions. 



2 Chemical Potential on the Lattice 

The chemical potential, /i, is introduced in the fermion action, i/jWiIj, as 0], 
W{x,x') = <5,_,,-«:^{(l-7,)[/,(x)5,,,^. + (l + 7,)[/t(x')5.,._j} 

i=l 

{e+^'^(l - 74)f/4(x)^.. .+4 + e-^'il + l4)Ulix')S.,,^,_i} . (1) 

Little is known about the behavior of dynamical fermion simulations when the chemical 
potential is introduced. For /i 7^ 0, the relation W'^ = 75^^75 does not hold, and hence 
detM^ is in general not real. 

Since U* = o"2?7^o"2, the fermion matrix for SU(2) has the following propertiy: 

W{x, x'; = - E {(1 - l*mx)5^,^,^i + (1 + i:)Uhx')6^,^,^.} 

1=1 

_^ {e+Ma(i _ 7*)[/4(x)5^, + e-'^^il + ll)Ul{x')5,,^,_i})a2, 
= a2W{x,x';j;)a2. (2) 

Then {detVF(x, x'; 'j^)}* = detVF(a:, x'; 7*). 7* belong a representation which also satisfies 
the anti-commutation relations same as 7^ and detVF should not depend on the represen- 
tation of 7-matrix. Therefore, differing essentially from the SU(3) case, the action of the 
SU(2)-QCD is real with chemical potential. 

However, numerical simulation is not straightforward and instability occurs with large 
chemical potential, which makes lattice simulation difficult Therefore, we need careful 
treatment for the updation of the configulation. We here adopt locally updating exact 
algorithm based on the Woodbery formula 0. The algorithm is summarized in the 
appendix. 



3 Thermo dynamical Quantities 

First, we calculate thermodynamical quantities, such as Polyakov line, gluon energy den- 
sity and baryon number density with 4^ x 12 lattice. We used Wilson fermion and Iwasaki 
improved action. 

The expectation value of the baryon number density is given by, 

<n>=-^^\ogZ (3) 

where Vg is the spatial volume Nx x Ny x Nz- Energy density, e, is given as. 
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Thermodynamical Quantities 




• Polyakov line 

♦ Gluon Energy Density *aM 
□ Number Density *a'^3 



Figure 1: Polyakov line, gluon energy density, and baryon number density as a function 
of chemical Potential //. 

The derivative of the partition function is composed of two parts, 

(logZ)' = i I VUV^l;Vi,{-S'a - S'p)e^-'^-'-\ (5) 

where Sq and Sp are gluon action and fermion action, respectively. We denote the 
contribution of the gluon action part (first term in the r.h.s. of eq.(3)) by gluon energy 
density. Figure 1 displays Polyakov line, gluon energy density and baryon number density 
of n where /3 = 0.7 and k — 0.15, respectively. All quantities start to have non-zero 
value at about fx = 0.4 and rise up with chemical potential in the region 0.4 < fj, < 0.8. 
None of them is order parameter of the phase transition in the exact sense, and since the 
lattice size is small, no sharp change is seen. However, growing up of the these quantities 
indicates that quarks and gluons become free from the confinement force at finite chemical 
potential about OA < fj, < 0.8. 

4 Meson and Baryon 

Because fundamental representation of the SU(2) is 2, baryon in the SU(2)-QCD is di- 
quark state. Scalar diquark state and pseudo-scalar diquark state are given as (C~^-?/')"^75-?/' 
and {C~^iIj)'^1'i/j, respectively, with C being charge conjugation matrix. Charge conjuga- 
tion makes transform property of the diquark state look opposite to the ordinary combi- 
nation of the scalar -iplip and pseudo scalar -ip^sip. We denote scalar diquark state and 
pseudo-scalar diquark state as b5b and bib, respectively, for the abbreviation. 
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Figure 2: Propagators of mesons and diquarks of 4^ x 12 lattice at /3 = 0.70 and k — 0.150. 
In order to avoid overplot of the data, diquark propagators are shifted 4 times. 



We evaluate propagators of the pseudo scalar iso vector meson, tt, the scalar iso vector 
meson, qq, vector meson, p, pseudo scalar diquark (bib) and scalar diquark (b5b). With 
vanishing chemical potential, the correlator of the tt and scalar diquark degenerate and 
so do Oq and pseudo scalar diquark. Because meson has no net baryon number, effect 
of the chemical potential is expected to appear in the meson propagator only through 
the mass. On the other hands, baryon (diquark) has definite baryon number; Hence, 
in addition to the change of the mass, the affect of the finite chemical potential on the 
particle and antiparticle is in opposite sign, which causes the asymmetry of the propagator 
in time direction. With a finite chemical potential /x, propagator of diquark, Gbu,{x,x') 



and Gb5b{x,x'), should behave as, 



Gb*b(t, /j) — Gb*b(T — t, — /i), 
in contrast to the meson propagator, 

with T being lattice size in the time direction. 

Figure ^ shows propagators of mesons and diquarks. At = 0, as expected, pseudo 
scalar vr and scalar diquark (b5b) and scalar oq and pseudo scalar diquark coincide, re- 
spectively (in Fig.^, to avoid the overplot, diquark propagators are shifted by factor 4). 
With finite chemical potential, asymmetry of the diquark propagators in time direction 
and anti-time direction become stronger with fi. On the other hand, meson propagators 
keeps symmetric in rit. Hence, in our calculation, effect of the finite chemical potential 
works appropriately on the hadron propagators. Hands P], reported that with finite 
chemical potential, propagator of the scalar diquark (b5b) and pseudo scalar diquark 
(bib) become parallel and they proposed the interpretation that both diquark states give 
the same spectra and the difference corresponds to the diquark condensation. However, at 
least the present statistics, our results do not give the the diquark propagators in parallel. 



5 CONCLUDING REMARKS 

We present numerical study of SU(2)-QCD with the chemical potential on lattice with 
Wilson fermions. Although the lattice is not large, behaviors of the thermodynamical 
quantities suggest that we are at around the confinement / deconfinement phase transition. 
Though the change of the propagator as a function of the chemical potential // is almost 
consistent with the results of 0, however, the spectrum of the diquarks seem to be 
different. But we need more statistics to conclude definite results. Estimation of the mass 
as a function of the chemical potential in the chiral limit is now in progress. 

In our calculation, numerical convergence becomes worse and worse with larger chem- 
ical potential and at ;U = 1.0 vanishing determinant makes simulation break-down. We 
are adopting the algorithm based on the exact calculation of the ratio of fermion deter- 
minant. Therefore, we can analyze the change of the distribution of eigenvalues of the 
fermion matrix with finite chemical potential 0. Investigation of the physical meaning 
of the numerical instability and distribution of the eigenvalue are also our next task. 
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7 Appendix 



We adopt an algorithm where the ratio of the determinant 

detW{U + AU) 



detW{U) 



det{I + W{Uy^AW) (6) 



is evaluated explicitly at each Metropolis update process, U U + AU, where AW = 
W{U + AU) — W{U) [|l^. An essential ingredient of the algorithm is Woodbury formula, 



{W + AW)-^ = W-^ - W-^AW{I + W-^AW)-^W-\ (7) 

Suppose we update link variables t/^(x)s only on a subset H of whole lattice. Though 
AW 7^ only on H, Woodbury formula still holds on H, and in this case, we can 
get the ratio of the fermion determinant as far as t/^(a;)s are locally updated only inside 
H. We take a 2^ hypercube as H. When we move to the next hypercube, {W~^)h^s are 
initialized by CG method. 

We employ an algorithm which takes into account the ratio of fermion determinant 
exactly, and has large Markov step, but we suffer from numerical instability at about 
/i = 1.0 
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